function [Xsource_1,Xsource_2,Xsource_3]=Unfolding(Xsource_p)
Xsource_1 = tenmat(Xsource_p,1); 
Xsource_2 = tenmat(Xsource_p,2); 
Xsource_3 = tenmat(Xsource_p,3); 
end